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ABSTRACT 



Type la supernovae (SNe la) are one of the major tools to 
determine the cosmological parameters. Utilizing them 
as distance indicators, it is possible to geometrically sur- 
vey the universe. To this end, the intrinsic scatter in the 
luminosities of these events needs to be calibrated using 
empirical relations between observables. A theoretical 
explanation for these relations is still lacking and can only 
be provided by a sound understanding of the mechanism 
of SNe la. Recently there has been significant progress 
in modeling SNe la. We report on numerical simulations 
of the explosion process, compare the results with obser- 
vations of nearby SNe la, and discuss current uncertain- 
ties of the models. The presented simulations shed some 
light on the origin of the diversity of SNe la. Such sim- 
ulations will pave the way towards an understanding of 
SN la diversities and correlations of their properties and 
ultimately provide a tool to validate the cosmological im- 
plications of SN la distance measurements. 

Key words: type la supernovae; numerical simulations; 
hydrodynamics; turbulence; radiative transport; cosmol- 
ogy- 



1. INTRODUCTION 



Almost 80 years after Einstein's attempts to derive a 
cosmological model from his theory of General Relativ- 
ity, subsequently introducing and discarding the cosmo- 
logical constant, two cosmol ogical surveys (R iess et aH 
119981 iPerlmutter etail 119991) of the universe based on 
distance measurements of type la supernovae (SNe la) 
revolutionized the cosmological paradigm. They found 
that the Universe is currently in a phase of accelerated 
expansion. This revived the discussion of the cosmolog- 
ical constant. Independent determinations of cosmologi- 
cal parameters based on anisotropies in the Cosmic Mi- 
crowave Background radiations and on large-scale galaxy 
surveys provided independent confirmation of the SN la 
measurements. 



Yet the question of the applicability of SNe la as distance 
indicators is still not satisfactory answered. SNe la are re- 
markably uniform events by astrophysical standards,but 
evidently no standard candles. Only the calibration of 
the distance measurements according to empirical corre- 
lations between observables provides the necessary ac- 
curacy for the determination of cosmological parameters. 
A firm theoretical reasoning of such correlations is, how- 
ever, still lacking. 

It is still an open question whether the accelerated expan- 
sion of the Universe is caused by a simple cosmological 
constant or some more complicated constituent. There- 
fore one of the most fascinating problems of cosmology 
is to determine the equation of state of this "Dark En- 
ergy" component. SNe la may be a suitable tool for this 
task, but the accuracy of measurements required here can 
only be reached if systematical errors that may arise from 
calibration of the distance measurements are under con- 
trol. A sound understanding of the physical mechanism 
of SNe la is inevitable in this project. 

This poses a great challenge to SNe la theory since it re- 
quires self-consistent parameter-free modeling so that the 
predictive power of the models is sufficient to tackle ques- 
tions of the origin of diversities and correlations. 

Consequently, as a first step, an astrophysical model has 
to be built and implemented into a numerical simula- 
tion code. As will be explained below, the nature of 
relevant physical effects requires three-dimensional sim- 
ulations. This approach is notably different from ear- 
lier one-dimensional SN la models which parametrized 
physical effects and therefore could be tuned to fit 
the observational data. Although this approach pro- 
vided valuable hints to underlying physical effects, it 
cannot be used to constrain the correlations. Three- 
dimensional models have been implemented over the past 



sumptions (e.g. Hillebrandt et al., 


2000; 


Reinecke et all 


2002b; Gamezo et al., 2003; Calder et al. 


2004- Most of 



them led to viable explosions and the task is now to com- 
pare the results of these simulations with observations of 
nearby SNe la. Once a model matching the observational 
requirements has been established, it can be used to ex- 
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plore the diversity arising from variations of physical pa- 
rameters as suggested by progenitor evolution and envi- 
ronmental effects. The results of such a survey should 
then be analyzed to identify correlations in the character- 
istics. This may eventually provide a theoretical basis for 
distance measurement calibrations and a handle on pos- 
sible systematical errors in those. 



2. ASTROPHYSICAL MODEL 



2.1. Progenitor evolution and ignition 

The cornerstones of the astrophysical model of SNe la 
are set by two fundamental characteristics of these events. 
Evidently, SNe la belong to the most energetic cosmic ex- 
plosions. For a short period of time they can outshine an 
entire galaxy consisting of tens of billions of stars. As- 
suming SNe la originate form single stellar objects, only 
their gravitational binding energy, released in a collapse 
towards a compact object, or its nuclear energy, released 
in explosive reactions, come into consideration as pos- 
sible energy sources. In the particular case of SNe la 
no compact object is found in the remnant excluding the 
first possibility. The fact that no hydrogen is found in the 
spectra of SNe la provides a strong hint that the object 
undergoing the nuclear explosion may be a white dwarf 
(WD) star consisting of carbon and oxygen. 

A simple WD, however, is an inert object and thus the 
only way to introduce the necessary dynamics into the 
system is to assume it to be part of a binary system. 
Here the WD accretes matter - in the favored scenario 
from a non-degenerate companion - until it approaches 
the Chandrasekhar mass (the maximum mass that is sup- 
ported against gravitational collapse by the pressure of 
the degenerate electrons). At this point, the density at 
the center of the WD increases rapidly so that fusion of 
carbon ignites. This gives rise to a stage of convective 
carbon burning lasting for several hundred years. This 
phase is terminated by one or more small spatial regions 
undergoing a thermonuclear runaway, marking the birth 
of a thermonuclear flame and the onset of the explosion. 

The limitation the fuel available in the explosion to the 
Chandrasekhar mass (~1.4M Q ) makes the described sce- 
nario particularly favorable since it provides a natural 
explanation for the striking uniformity of SNe la in the 
gross observational features. On the other hand, it is 
afflicted with great uncertainties. Achieving a stable 
mass transfer in the progenitor binary system to build 
up a Ch andrasekhar mass WD is highly non-trivial (e.g. 

and the observational evidence 
for such systems is sparse. The convective burning stage 
and the conditions at flame ignition are extremely hard to 
model both analytically and numerically. Therefore the 
exact shape and configuration of t he first flame spark(s) 
is not yet well const rained (e.g. IWooslev et all 12004 
iHoflich A Stei r H200l . 



2.2. Flame propagation and explosion 

The goal of SN la explosion models is to follow the prop- 
agation of the thermonuclear flame from its ignition near 
the center of the WD outwards and to determine the com- 
position and the distribution of the burning products in 
the ejected material. 

In principle, hydrodynamics allows for two different 
modes of flame propagation. One is the subsonic de- 
flagration in which the flame is mediated by the thermal 
conduction of the degenerate electron gas and the other 
is a supersonic detonation in which the flame is driven 
by shoc k waves. Alr eady the first attempts to simulate 
SNe la ( lArnettt 1 19691) revealed that a prompt detonation 
cannot explain these events since the full star would be 
incinerated with sound speed. All material would then 
be burnt at high densitites and converted to iron group 
elements. Observed spectra of SNe la, however, clearly 
show indication for significant amounts of intermediate 
mass elements (like Si, Ca, S . . . ) in the ejecta. These 
can only be produced in the burning at lower densities 
and therefore the material needs to expand prior to incin- 
eration. This is provided if the flame propagation starts 
out as a slow deflagration. A laminar deflagration front, 
however, propagates too slowly to release sufficient en- 
ergy to explode the WD star. 

Thus, a key feature of a viable SN la model is the way of 
flame acceleration. For this, two effects are concievable. 

Firstly, the flame propagation from the center of the 
star outwards is buoyancy unstable, since it leaves be- 
hind light and hot ashes below the dense fuel - a den- 
sity statification inverse to the gravitational acceleration. 
In its non-linear stage, the Rayleigh-Taylor instability 
leads to the formation of mushroom-shaped burning bub- 
bles raising into the fuel. The Reynolds number typi- 
cal for this situation is as high as 10 14 . Clearly, shear 
(Kelvin-Helmholtz) instabilities at the interfaces of these 
bubbles will generate turbulent eddies which then de- 
cay to smaller scales foming a turbulent energy cascade. 
The flame will interact with these eddies down to the 
Gibson-scale at which the turbulent velocity fluctuations 
become comparable to the laminar flame speed. This in- 
teraction corrugates the flame increasing its surface and 
consequently enhances the net burning rate. Thus, the 
flame is significantly accelerated and three-dimensional 
simulations (e.^. fHillebran dt et allfeoOOtlReinecke et all 
2002b: lGamezo et alll200 3) show that this effect may in- 
deed lead to an explosion of the WD. In this deflagration 
model, SN la explosions are a problem of turbulent com- 
bustion. 

Secondly, a transition of the flame propagation mode 
from subsonic deflagration to supersonic detonation 
could provide an additional enhancement of the burning. 
However, a mechanism providing such a transition in a 
SN la explosion could not be found yet and is therefore 
hypothetical. It can only be introduced at hoc into a SN 
la model and is therefore free parameter - an effect that 
one would rather avoid in order to maintain the predictive 
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power of the model. 

We will thus focus on the pure deflagration model in the 
following. After outlining the numerical implementation 
in Sect. [3] we will thus discuss the results from simula- 
tions and the question to which degree they meet observ- 
able constraints in Sect. |4] 



3. NUMERICAL IMPLEMENTATION 



3.1. Relevant scales 

The numerical implementation of a deflagration SN la 
model is significantly complicated by the wide range of 
relevant length scales involved in the problem. From the 
radius of the WD star (~2000km at the onset of the ex- 
plosion and expanding in the process) it reaches down to 
the flame width and the Kolmogorov scale at which the 
turbulent energy is dissipated - both well below one cen- 
timeter. Simulations capturing the entire star reach reso- 
lutions around one kilometer while the Gibson scale is of 
the order of 10 4 cm at the beginning of the explosion and 
decreases steadily. 

For large-scale SN la simulations this has three conse- 
quences. 

1 . It is not possible to fully resolve the interaction of 
the flame with turbulence. Therefore modeling of 
the effects on unresolved scales is necessary. 

2. The internal flame structure cannot be resolved. 
Compleme ntary sma l l scal e simulatio ns in one 
llTimmes & Wooslevl 1 19921) or more JBell et all 
2004) dimensions are required here. 

3. Assumptions about the flame properties at unre- 
solved scales (e.g. stability below the Gibson scale) 
have to be validate d in separate sm all-scale simula- 
tions (iRopke et alll2003ll2004albl) . 



3.2. Numerical techniques 

As obvious from the consideration of the relevant scales, 
a suitable approach to implement a SNe la deflagra- 
tion model on scales of the WD is that of a Large 
Eddy Simulation (LES). It captures the effects of tur- 
bulent phenomena at the largest scales where energy is 
fed into the turbulent cascade and models turbulence 
on unresolved scales. The numerical techniques form- 
ing the foundation of the implementation of a defla- 
gration SN la model d e scribed here a r e described in 
iReinecke et all lll999bl). iRopkd J2005I). ISchmidt et all 
l l2005l) . and lNiemever & Hillebrandtl i ll 9951) . 

The hydrodynamics on the resolved scales is modeled in 
a finite volume approach. It is based on the Prometheus 



implementation dFrvxell et all 119891) of the piecewise 
parabolic method (Colella & Woodward, 1984). 

As a consequence of \l\ of the enumeration above, 
a sub grid- scale model is applied to account for tur- 
bulent effects on unresolved scales. Some older 
simulation s follow the implementation suggested by 
iNiemever & Hillebrandtl dl995l) . while one recent highly 
resolved run applied the updated modeling approach of 
ISchmidt et all J2005h . 

According to[2]from above, the flame representation has 
to follow a modeling approach. Seen from the scales 
of the WD star, it appears as a sharp discontinuity sep- 
arating the fuel from the ashes. A suitable numerical 
method to follow the evolution of such an interface is 
the level set tec hnique introduced to SN la modeling by 
IReinecke et ail jl999al) . 

In such an implementation, the flame propagation speed 
needs to be prescribed. This is, however, not an arbi- 
trary quantity in our modeling framework. For turbulent 
combustion in the flamelet regime, which applies to the 
burning in major parts of SN la explosions, it is given by 
the turbulent velocity fluctuations. These are determined 
by the subgrid-scale turbulence model. 

Nuclear reactions are implemented i n the si mplified ap- 
proach proposed by Reinecke et all y002a). The pro- 
genitor material is composed of 12 C and 16 0. At high 
fuel densities nuclear burning terminates in nuclear sta- 
tistical equilibrium represented by a mixture of 56 Ni and 
a-particles. Once the fuel density drops below 5.25 x 
10 7 g cirT 3 burning will become incomplete and interme- 
diate mass elements (represented by 24 Mg) are produced. 
The respective differences in nuclear binding energy is 
released which provides sufficient accuracy to model the 
dynamics of the explosion. 

However, in order to derive synthetic spectra and light 
curves from models, the accurate distribution of different 
chemical species in the ejecta needs to be determined. 
This is achieved in a postprocessing step which makes 
use of the temperature and density evolution recorded by 
tracer particles advected in the hydrodynamical explosion 
simulation. On the basis of this data, the details of the nu- 
clear reactions in the ashes are reconstr ucted emp loying 
an ex tended nuclear reaction network ( Travagl io et all 
2004). 



4. SIMULATIONS VS. OBSERVATIONS 



4.1. Explosion models and global characteristics 

We will now review the status of deflagration SN la mod- 
eling with particular emphasis on the question how the 
results of such simulations compare with observations. 
Due to recent progress in such simulations, a direct com- 
parison with details of observations of nearby SNe la has 
come into reach. 



Figure 1. Snapshots from a high-resolved SN la simulation starting from a multi-spot ignition scenario. The logarithm of 
the density is volume rendered indicating the extend of the WD star and the isosurface corresponds to the thermonuclear 
flame. 



One requirement to reliably derive observables from sim- 
ulations is that the evolution of the models be followed to 
the stage of homologous expansion. In this hydrodynam- 
ically relaxed situation, the velocity of the ejected mate- 
rial is proportional to its radius. The implementation of a 
moving computational grid on which the simulations are 
carried out facilita ted the long-t erm evolution of expand- 
ing SN la models (RopkeJUOQH). 

Due to the high computational expenses, most three- 
dimensional simulations comprised only one spatial oc- 
tant of the WD star, assuming mirror symmetry to the 
other octant s. However, only full-star setups (such as 
presented by Ro pke & Hillebrandtl (120051) 1 allow to study 
asymmetry effects. On the basis of spectrapolarimetry 
observations of several SNe la (e.g. IWang et all 120031) 
these are expected to occur in at least some explosions. 
Ropke & Hillebrandt ( 2005 ) showed, that such asymme- 
tries arise exclusively from asymmetries in the flame igni- 
tion conditions and not from large-scale instabilities and 
resulting prefered modes in the flow patterns. Therefore, 
fixing a symmetric initial flame shape and studying the 
influence of other physical parameters on single-octant 
explosion models is a valid approach. 

The propagation of the thermonuclear flame through the 
WD is shown in Fig.^ The snapshots are taken from a 
recent high-resolution simulation in which the flame was 
ignited in multiple sparks around the center of the star (cf. 
the left part of Fig. It agrees with the expectations in 
the astrophysical scenario outlined in Sect. 12.21 After ig- 
nition, the Rayleigh-Taylor instability leads to the forma- 
tion of mushroom-shaped burning bubbles on large scales 
which are subsequently superposed by smaller modes of 
the instability. The resulting flow quickly becomes turbu- 
lent and the flame is being corrugated. This is apparent 
in the second snapshot of Fig.^ In this turbulent burning 
mode, the flame propagates through most parts of the star 
(cf. the third snapshot in Fig.^. A powerful explosion of 
the WD star results in this simulation. 



Naturally, the question arises if the outcome of simu- 
lations as the one described above meets the observa- 
tional constraints. Such constraints arise from the global 
characteristics derived from observations, observed light 
curves, and spectra taken from nearby SNe la. 

The global characteristics derived from SN la observa- 
tions state that a valid explosion model should release 
around 10 51 erg of energy and produce ~0.4...0.7M o 
of 56 Ni in the nuclear burning. However, there exists 
a large diversity in the observations ranging from the 
class of sub-luminous SNe la (like SN 1991bg) to super- 
luminous events (e.g. SN 1991T). Our deflagration mod- 
els started with a multi-spot ignition setup typically pos- 
sess 7 ... 8 x 10 50 erg of asymptotic kinetic energy of the 
ejecta. They produce around 0.4 M of 56 Ni. Thus they 
fall into the range of observational expectations, but in 
the current stage do not account for the more energetic 
SNe la. 



4.2. Lightcurves 

A further requirement is that synthetic lightcurves agree 
with observed ones. These are sensitive to the energy 
release, the 56 Ni production, as well as to the distribu- 
tion of elements in the ejecta. In Fig. [2] we compare 
synt hetic light curves derived from the 2_2_2 simula- 
tion fitopke et ali 120051) with observations of SN 1994D. 
The multi-band light curve of this model was ca lcu- 
lated using the Stella c ode of Blin nikov et alJ Jl99 8 ) and 
iBlinnikov & Sorokinal ( l2000l) . 

The distance to SN 1994D is still controversial. 
We have assumed the dist ance modulus 30.4 from 
iDrenkhahn & RichtleJ(ll999l) 

The model produced 0.3 M of 56 Ni, although observa- 
tions require a bit higher 56 Ni mass (~ 0.4 M ) for the 
assumed distance. There is generally very good agree- 
ment in the B and V bands near peak luminosities and in 
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Fi gure 2. Synthetic lig ht curves derived from model 2_2_2 
of (Ro oke et all 1200.51) (solid curves) compared with ob- 
served light curves from SN 1994D. 



decline rate 20 days after the peak which is most impor- 
tant for cosmological applications of type la supernovae. 



4.3. Spectra 

A much harder test for the models is posed by the com- 
parison of synthetic with observed spectra since these de- 
pend on details in the compositon of the ejected material. 

Nebular spectra provide a means of studying the central 
parts of the ejecta, since they are taken at epochs where 
these have become transparent due to expansion. Thus 
they explore the "heart of the explosion" and are a valu- 
able tool to study details of the physical processes in- 
volved in the explosion stage. Unfortunately, only one 
single synthetic late time sp ectrum is ava i lable from our 
deflagration SN la models ( Kozm a et all 120051) . It was 
derived from a run that was intended to test the imple- 
mentation of the moving grid and is not believed to pro- 
vide a realistic SN la model. The flame was ignited in 
an axisymmetric shape where a sphere around the center 
of the WD star was perturbed by a toroidal structure of 
large wavelength and amplitude. This artificial and sim- 
ple initial flame shape gives reason to not expect a good 
agreement between model and observation. However, as 
shown in Fig. [3] the broad iron features of the observed 
spectra are qualitatively reproduced. An inconsistency of 
the model with the observed nebular spectra is the ap- 
pearance of a pronounced oxygen line at 6300 A. Both 
features of the synthetic spectrum share a common ori- 
gin. The broad iron lines are caused by NSE material that 
is transported in the uprising plumes of ashes and thus 
distributed in velocity space. At the same time strong 
downdrafts carry unburnt material towards the center of 



The disagreement may in part be attributed to the sim- 
plicity of the explosion model. Its highly symmetric ini- 
tial flame shape with large imprinted perturbations favor a 
pronouned evolution of large-scale Rayleigh-Taylor fea- 
tures. To answer the question whether the unburnt mate- 
rial at low velocities is a generic problem in the deflagra- 
tion model of SNe la, the high-resolved simulations men- 
tioned above was carried out. The finer structure of the 
turbulent flame in this simulation together with improved 
flame modeling led to a better burnout of the downdrafts. 
Much less unburnt material is found at low velocities in 
this model. 



A powerful diagnostic tool to compare SN la models with 
observations is provided by th e abundance tomography 
presented by Stehle et al. (2005). It makes use of spectra 
taken from SN 2002bo with an extraordinary good time 
coverage. Fitting this sequence of data with synthetic 
spectra unveils the composition of the ejecta in velocity 
space slice by slice, since the photospere moves gradu- 
ally inwards with the expansion of the remnant. This 
abundance tomography of the ejecta can be compared 
with results of our three-dimensional models, when av- 
eraged over the angles. Qual itatively, the mixed com- 
position of the ejecta found bv lStehle et alJ ( 120051) is re- 
produced by deflagration SN la models in a natural way 
since these predict a distribution of burnt material with 
the rising bubbles. A problem was, however, that older 
predicted large unburnt material fractions in the central 
parts of the ejecta in disagreement with the results of 
Stehle et al. (2005). Our recent high-resolved simulation 
cures this problem by clearly reproducing the iron-group 
dominance in the low-velocity ejecta. The results are cur- 
rently being analyzed and seem promising for matching 
the observations. 
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Table 1. Variation of initial parameters in SN la explosion 
models. 



Parameter 


range of 
variation 


effect on 
56 Ni pro- 
duction 


effect on 
total en- 
ergy 


X( 12 C) 


[0.30,0.62] 


<2% 


-14% 


Pc [10 9 g/cm 3 ] 


[1.0,2.6] 


-6% 


-17% 


Z[Z Q ] 


[0.5,3.0] 


-20% 


none 



5. DIVERSITY AND CORRELATIONS 



The recent developments in the deflagration SN la ex- 
plosion modeling outlined in the previous section seem 
to indicate that such a model is capable of reproducing 
the main features of observed objects. They do not rule 
out the alternative of a delayed detonation. The success 
of the pure deflagration model, however, may be taken 
as evidence that this stage provides at least a major con- 
tribution to the explosion process and general properties 
will not be changed much by a hypothetical deflagration- 
to-detonation transition. 

It is therefore an important question to ask how such a 
model behaves under variation of physical parameters. 
Does it reproduce the observed diversity of SNe la? Are 
correlations between observables evident in the results? 

Unfortunately, three-dimensional deflagration models of 
SNe la as described above are computationally expen- 
sive. To moderate these expenses, simplified setups may 
be used to study effects of physical parameters on the ex- 
pl osion models . Such an approach was recently taken 
by Ropke et al. (2005) and resulted in the first system- 
atic study of progenitor parameters in three-dimensional 
models. The basis of this study was a single-octant setup 
with moderate (yet numerically converged) resolution. 
However, the lack in resolution did not allow a reasonable 
multi-spot ignition scenario and thus only weak explo- 
sions can be expected. It is therefore not possible to set 
the absolute scale of effects in this approach, but trends 
can clearly be identified. 

The parameters chosen for the study were the WD's 
carbon-to-oxygen ratio, its central density at ignition and 
its 22 Ne mass fraction resulting from the metallicity of the 
progenitor. All parameters were varied independently to 
study the individual effects on the explosion process. In 
a realistic scenario, however, these parameters are inter- 
related by the evolution of the progenitor binary system. 
The results of this survey are given in Tab. [2 

A variation of the carbon-to-oxygen ratio effects the en- 
ergy production in the burning due to the differences in 
the binding energies of these two nuclei. Counterintu- 
itively, this results in no significant cha nge in the 56 Ni 
production. Ropke & Hillebrandt (2004) explained this 



effect by the fact that the potentially higher energy release 
in carbon-rich models is buffered by a higher a-particle 
fraction in the NSE and only released when burning is 
already incomplete. 

Models with lower central densities show a delayed flame 
evolution and consequently a lower and delayed energy 
production. This is due to the fact that the flame expe- 
riences a lower gravitational acceleration in these mod- 
els resulting in decreased turbulence generation. There- 
fore less 56 Ni is produced in these models. This effect 
is even more pronounced due to the fact that in low- 
density WDs less material is present at sufficient densi- 
ties to be potentially burnt to NSE. A counteractive ef- 
fect is expected at higher densities. Here, electron cap- 
tures will become important favoring neutron-rich iso- 
topes over 56 Ni in the NSE. The dynamical effects of 
electron captures are, however, not yet implemented in 
our explosion model and therefore the current survey is 
restricted to relatively low central densities. 

The metallicity of the progenitor star results in a certain 
22 Ne mass fraction in the WD. This is an isotope with 
neutron excess and therefore again favors the production 
of neutron-rich species over 56 Ni in the NSE. Our results 
confirm the an alytic prediction bv iTimmes et alJ ll2003l) 
and agree with Travariio et aD ( 120051) . The metallicity 
parameter, however, has no effect on the explosion dy- 
namics and the energy production in our models. 

To determine the effects of these variations on observ- 
ables, synthetic light curves were derived from all mod- 
els (an example is shown in Fig. [5). From these, the peak 
luminosities and decline rates (in magnitudes 15 days af- 
ter maxi mum; Amis), w ere determined. The pioneering 
work by Phillips! dl993l) established the relation between 
the peak luminosity and Amis as one of the primary tools 
to calibrate cosmological SN la distance measurements. 
The so-called Phillips relation quantifies the decrease of 
Amis for brighter SNe la. 

The results fr om our model s are co mpared with the rela- 
tion given by Phillips et al. ( 1999) in Fig.0] Obviously, 
the absolute magnitude of the Phillips et al. ( 1999) rela- 
tion is not met by our set of models (cf. the upper panel 
of Fig. |4}. Moreover, the range of scatter in Amis is 
m uch narrower than that of the set of observations used 
by Phillips et al] d!999h . The trend of our models is not 
clear. T hey are best fit by a line with a slope consistent 
with the Phillips et al J 1119 99) relation. It is, however, ob- 
vious, that a better agreement cannot be expected given 
the fact, that the parameters in our survey were chosen in- 
dependently. A consistent stellar evolution would pick a 
sub-sample of our set of models and it will be investigated 
in forthcoming studies whether this will be in agreement 
with the Phillips relation. 

Nonetheless, with our set of models, the question can be 
answered, which parameter dominates the slope in the di- 
rection of dimmer events for faster decline rates. The var- 
ied parameters are coded by different line styles in the 
lower panel of Fig.|4] Clearly, the progenitor's metallic- 



7 



o« o. 



1.4 1.5 1.6 1.7 

Am 15 (S) 



-18.3 - 



=1 



-18.2 - 



-18.1 - 




Figure 4. Peak luminosity vs. decline rate of the light 
curve in the B band ( diamonds correspond to SN la ex- 
plosion models). C ompared with original relation by 
Philli ps "et al\ fc^99) (dashed curve) and shifted relation 
(solid curve) in the top panel. 



ity can be identified as this parameter. Variations in the 
central density and the carbon mass fraction of the WD 
superpose a scatter on the dominant relation. 



6. CONCLUSIONS 



We discussed the mechanism of SN la explosions and 
reviewed the common astrophysical scenarios applied to 
explain these events. The implementation of such models 
into numerical simulations was described focusing on the 
pure deflagration model for SNe la. Subsequently, recent 
progress in numerical simulations was reported. On the 



basis of comparison of these results with observations, 
the conclusion is drawn, that pure deflagration models 
are at least capable of reproducing the gross features of 
observed SNe la. 

Such models therefore can be applied to study the origin 
of the diversity of SNe la. In a first systematic survey 
of initial parameters of three-dimensional models, the ef- 
fects of the WD's carbon-to-oxygen ratio, its central den- 
sity and the progenitor's metallicity on the explosion pro- 
cess were analyzed. Although this was only possible on 
the basis of simplified and therefore weakly exploding 
models, clear trends could be identified. 

From this set of models, synthetic light curves were de- 
rived. The relation between the peak-luminosities and the 
decline rates of the light curves were found to be con - 
sistent with the slope predicted by Phillies et al. ( 1999). 
Due to the simplicity of the applied models, the absolute 
magnitude of effects could not be predicted. 

Forthcoming studies will focus on the questions how a 
realistic stellar evolution of the progenitor interrelates 
the effects. A second important question is whether the 
range of scatter predicted by the observations can be re- 
produced with more elaborate models. Although it cap- 
tures main aspects of observations, it may be possible 
that the applied SN la model is still incomplete and addi- 
tional physical effects (like delayed detonations) need to 
be taken into account. 
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